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We compute the one-loop self-energy correction to the massive domain wall quark propagator. Combining this 
calculation with simulations at several gauge couplings, we estimate the strange quark mass in the continuum 
limit. The perturbative one-loop mass renormalization is comparable to that for Wilson quarks and considerably 
smaller than that for Kogut-Susskind quarks. Also, scaling violations appear mild in comparison to other errors 
at present. Given their good chiral behavior and these features, domain wall quarks are attractive for evaluating 
the light quark masses. Our preliminary quenched result is ms{2 GeV) = 82(15) MeV in the MS scheme. 



1. INTRODUCTION 

Computing light quark masses is a high prior- 
ity. Present lattice predictions of mi = (m„ + 
md)/2 and rris give differing results depending 
on the particular method. Domain wall (DW) 
fermions ||l|,|| respect the chiral symmetries of the 
continuum exactly in the limit Ns —>■ oo, be- 
ing the number of sites in a fictitious extra dimen- 
sion. This is an especially attractive feature for 
simulating light quark physics where chiral sym- 
metry is crucial 1^,^. Furthermore, simulations 
have demonstrated that Ng ^ 10 for /? > 6.0 
suffices for very good chiral behavior rendering 
DW fermions quite practical Q. We report on 
first results for the strange quark mass using DW 
fermions. 

2. ACTION 

We use the boundary fermion variant of the 
DW formulation. The fifth dimension has a finite 
number of sites ; a light right-handed mode is 
bound to the 4-d surface at one end and a light 
left-handed mode at the other end. A 5-d mass 
M determines the strength of the binding, and a 
parameter m controls the coupling between the 
two 5-d boundaries. 

The chiral symmetry breaking is due to the 
explicit term proportional to m and to implicit 
mixing between the two modes which should be 
suppressed ~ e~"^' , where a > 0. Therefore the 



latter can be made arbitrarily small compared to 
the former by increasing Ns- For example, at tree 
level the quark mass is 
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If, in the free theory, M is in the range < 
M < 2, there is a single light mode fixed to either 
boundary. More properties of DW fermions are 
discussed in |^ and references therein. 

3. RENORMALIZATION 

Interactions renormalize the five-dimensional 
mass M additively; thus the range < Af < 2 is 
shifted by an amount Mc. In perturbation theory 
Mc can be computed from terms proportional to 
a~^. Numerically, the tadpole graph gives larger 
contributions to the self-energy than the half- 
circle graph; so one way to estimate Mc is to use 
the tadpole term alone, giving Af''"^ — 12. Ga^ 0- 
Noting the importance of Mc and the limitations 
of perturbation theory, a nonperturbative deter- 
mination of Mc is more desirable. Through the 
overlap formulation |^ the Wilson-Dirac opera- 
tor defines a transfer matrix which governs prop- 
agation in the 5th dimension. Therefore Mc can 
be determined directly from the critical Wilson 
hopping parameter M: 
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Using data from ||l^ we give M^ relevant to this 
work in Table |^. The superscript is to distinguish 
this calculation of Mc from other estimates. For 
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Figure 1. The matching coefficient Cm- The 
dashed lines indicate the range of M — Mc where 
simulations have been performed. 



Table 1 

Summary of simulation parameters and results. 



,LAT 



and mf ^(2 GeV) are in MeV. 





/? = 6.3 


(3 = 6.0 


(3 = 5.85 


# configs. 


18 


33 


18 


volume 


243 X 60 


16^x32 


16^x32 




10 


10 


14 


M 


1.5 


1.7 


1.7 




0.708 


0.819 


0.908 


Cm 


2.35 


2.47 


2.42 


(GeV) 


3.4(3) 


2.39(12) 


1.69(10) 


^1 strange 


0.020(3) 


0.027(3) 


0.041(6) 


777,LAT 


64(11) 


62(8) 


66(10) 


(Tr C/piaq/3) 


0.6224 


0.5937 


0.5751 




0.131 


0.146 


0.157 


S(2 GeV) 


85(15) 


81(10) 


83(13) 



reasonable definitions of the coupling constant as 
(see later), we find M^''^'^ « . 

In order to compute the quark mass renormal- 
ization, we have extended the explicit one-loop 
calculation of the DW fermion self-energy to 
nonzero fermion mass jll]] . The renormalization 
factor Zm can be defined by equating the one- 
loop fermion propagators on the lattice and the 
continuum |l|]. The result IS 



Z„i = 1 — {2as/TT) ln(/ia) — C, 



(3) 



where the DW quark masses are computed with 
a lattice spacing a and the MS quark masses at a 
scale /i. Cm depends only on M — Mc and is plot- 
ted in Fig. |l|. The Cm for the M - used in 
our simulations are shown in Table ^ and should 
be compared to 2.16, 3.22, and 6.54 for Wil- 
son, Sheikholeslami-Wohlert (SW), and Kogut- 
Susskind (KS) fermions. 

4. QUARK MASS 

We have computed the pion mass and decay 
constant on a few dozen quenched configurations 
at /? = 5.85, 6.0, and 6.3. We give other simula- 
tion parameters in Table 0. In addition we stud- 
ied the pion mass as a function of both m and M 
a,t (3 — 6.0. For each M, the pion mass squared 



extrapolates linearly to zero at m = within er- 
rors as in ^ . Fig. || shows the dependence of the 
pion mass squared on M. In lowest order chiral 
perturbation theory ~ niq. Therefore, from 
Eqn. (0) with M ^ M - the pion mass 

squared should obey 

(aM,)2 = (const.) (M - M^){2 - M + M]^).(4) 



The dotted lines are fits to (Q) and have good 
X^'s. At present the data do not permit a more 
general fit. 

To numerically compute the renormalization 
constant (^) , one must choose a definition for the 
coupling constant as and the scale ^ at which 
it is evaluated. We follow ||l^ and use a-^{iJ,). 
Specifically, we compute ay (3.41 /a) from the pla- 
quette and convert to perturbatively |l^ 

We run that coupling constant to to apply (| 
and finally run m-^ perturbatively to 2 GeV. We 
choose fj, — 2/ a for the matching; the final re- 
sult varies by 2 MeV for 0.5 < /ia < tt. We use 
the physical pion decay constant to set the lat- 
tice spacing and the physical kaon mass to fix the 
parameter m to the strange quark system. 

In Fig. ^ we plot our final results for the strange 
quark mass along with recent results using other 
actions. The results for DW fermions appear 
to be in rough agreement with other lattice cal- 
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Figure 2. Pion mass squared as a function of 
M at /? = 6.0. N, = 14 except for the M = 1.7, 
m = 0.075 point where Ng = 10. The dotted hues 
are fits to Eqn. H. 
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Figure 3. Strange quark mass in MS scheme 
at 2 GeV. Our results are displayed as squares, 
Wilson fermions as diamonds ||l5|, SW fermions 
as circles Ml, and KS results as asterisks |17|,|l8|. 



culations of the strange quark mass using the 
kaon mass and the perturbative Within our 

rather large errors, the DW results appear scale 
independent, so for now we take a weighted av- 
erage. Our statistical errors are large due to our 
small data set. We attribute a 2 MeV error to the 
ambiguity in q* and tentatively assign a 15% error 
due to determining from /tt Quenching 
errors are not included at present. We add the 
systematic and statistical errors in quadrature to 
obtain the preliminary result for the MS quark 
mass: to, (2 GeV) = 82(15) MeV. 

We find that DW fermions are a viable method 
of studying the light quark masses. We have 
shown that perturbative corrections to the mass 
renormalization are under control, the lattice 
quark mass depends on to and M as expected, 
and scaling violations appear mild at this stage. 

We thank T. Bhattacharya for discussions. 
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